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I 


Introduction 


An  electromagnetic  wave  propagating  through  the  earth's 
atmosphere  is  subject  to  attenuation  due  to  scattering  by  molecules, 
aerosols  and  turbulence  (which  acts  like  a  source  of  inhomogeneities) 
as  well  as  due  to  absorption  of  radiation  by  atmospheric  constituents 
in  the  atmosphere.  In  addition  to  attenuation,  scattering  by 
turbulence,  in  which  the  refractive  index  varies  through  space  and 
time,  makes  the  wave  intensity  fluctuate,  especially  for  waves  with 
wavelength  shorter  than  millimeters.  Because  the  wavefront  is  spread 
by  the  smaller  scales  and  wandered  by  the  larger  scales  of  the 
turbulence,  the  reformed  wavefronts  cause  the  scintillation  of  the 
wave  intensity. 

In  this  report,  we  review  theoretically  the  attenuation  and 
fluctuation  of  the  wave  propagating  through  turbulent  atmosphere  and 
Incorporate  the  results  into  the  Lowtran1,  computer  code  in  terms  of 
two  new  subroutines.  These  added  subroutines  are  used  to  calculate 
the  normalized  intensity  (or  power)  variance,  by  which  we  define  the 
upper  and  lower  bounds  of  atmospheric  transmittance.  Since  the  wave 
attenuation  due  to  turbulence,  which  will  be  discussed  in  section  II, 
is  not  significant,  we  do  not  Incorporate  it  into  the  calculation  of 
atmospheric  transmittance  in  Lowtran  in  which  the  attenuation  due  to 
molecular  and  aerosol's  absorption  and  scattering  is  taken  into 
account. 

In  practical  systems,  transmitters  and  receivers  with  finite 
apertures  are  used.  A  larger  receiver  aperture  not  only  collects  more 
power,  but  also  reduces  wave  fluctuation.  This  is  called  the  receiver 
aperture  averaging  effect.  Also,  a  more  coherent  source  gives  rise  to 


smaller  scintillations.  Hence,  the  effects  of  aperture  size,  source 
coherence  and  turbulence  on  wave  scintillation  are  included  in  the 
formulation  which  will  be  used  to  code  the  new  subroutines  in 


Lowtran. 


We  have  chosen  pair-correlated  field  statistics  ,  which  act  like 
Gaussian  field  statistics,  to  model  the  coherence  poperties  of 
sources.  It  is  shown  that  these  statistics  yield  satisfactory 

results^. 

The  extended  Huygens-Fresnel  method"*  has  been  used  to  obtain 
scintillation  expressions  for  partially  coherent  beam  waves  as  well  as 
for  spherical  waves.  Since  the  wave  structure  functions  which  are 
used  are  valid  for  aperture  sizes  which  are  smaller  than  the  Fresnel 
zone  and  since  the  parallel  approximation  is  applied  in  the 
Huygens-Fresnel  formulation,  we  cannot  extend  the  beam  wave  result  to 
the  plane  wave  case  simply  by  letting  transmitter  size  go  to  infinity. 
Therefore,  we  derive  the  intensity  variance  for  plane  wave  sources  by 
use  of  the  log-amplitude  variance  which  can  be  obtained  by  Rytov's 

method**.  Thus,  two  subroutines,  VRANI  and  SINTL,  are  coded  separately 
for  plane  waves  and  beam  waves  (including  spherical  waves), 
respectively. 


v-  *VV*  ■V*.  -  WvV 

V*  •  «  *a>_a  '•*  »>a 


II.  Transmittance  attenuation  due  to  turbulence 


The  transmittance  for  a  wave  propagating  through  the  atmosphere 
is  defined  as 


T 


<I> 


<IV> 


(1) 


v 

where  I  and  I  are  the  received  wave  intensities  in  the  atmosphere 
and  in  vacuum,  respectively.  <  >  denotes  an  ensemble-average. 
Excluding  the  attenuation  due  to  turbulence  scattering,  the 
atmospheric  transmittance  described  in  LOWTRAN  is 


TL  Tkm  *  Tka  *  T«n  *  TOa 


(2) 


where  Kn>  ,  a*  are  the  absorption  constant  and  scattering  cross 

section  of  the  molecules  and  aerosols  in  the  atmosphere,  respectively. 
Eq.(2)  implies  that  all  scattering  energy  is  thought  as  a  loss.  This 
is  a  good  approximation  for  a  receiver  with  a  very  narrow  angle  of 
Field  of  View  (FOV)  and  for  scattering  by  molecules  and  aerosols  that 
give  trivial  forward  scattering  because  of  their  small  sizes. 

However,  the  scales  of  turbulence  ( 1mm  ~  100m)  are  much  larger  than 
optical  or  infrared  wavelengths.  A  strong  forward  scattering  field 
due  to  turbulence  is  then  present.  Since  the  turbulence-induced 
deviation  of  the  arrival  angle  is  in  the  order  of  micrc -radians  ,  the 
scattering  energy  is  not  completely  lost  for  a  receiver  with  an  angle 
of  FOV  larger  than  micro- radians.  Assuming  that  the  scattering  by 
molecules,  aerosols  and  by  turbulence  are  Independent,  we  can  then 
write  the  atmospheric  transmittance  t  as 


flUM'J1 


»'  iR  jr  .  v*t  m.' 
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T  -  T  <T  > 
L  T 


(3) 


where  <r  >  is  the  ensemble-averaged  transmittance  due  to  turbulence. 

(a)  Plane  wave 

In  turbulence,  the  refractive  index  variation  is  small  (~  10*^) 
and  thep  scales  of  turbulence  are  much  larger  than  optical  or  infrared 
wavelengths.  Therefore,  the  backsca ttering  due  to  turbulence  is 
small.  For  plane  waves,  neglecting  the  trivial  backscattered  fields, 
almost  all  incident  waves  reach  the  receiver  plane  though  the 
wavefronts  are  distorted  and  the  received  wave  energy  is  redistributed 
through  the  entire  receiver  plane.  Hence,  the  ensemble -averaged 
transmittance  of  a  plane  wave  in  turbulence  has  a  value  of  unity, 
i.e. , 


<Y 


1  (4) 

However,  the  redistributed  energy  does  degrade  the  coherence  of  the 
received  wave  and  induces  wave  scintillation  which  we  will  discuss  in 
section  III. 

The  finite  aperture  of  a  receiver  with  radius  R  collects  power  P: 


<P>  -  I  tuR 
s 


(5a) 


■  t,I  tR 
L  s 


(5b) 


where  1^  is  the  plane  wave  intensity  in  the  transmitter  plane.  The 
transmittance  of  power  can  then  be  defined  as 


<P>_ 

"*v" 

<P  > 


T  <  T  > 

L  p,T 


(6) 
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It  Is  obvious  the  thet  <PV>  -  I  *R2.  Substituting  <PV>  into  Eq.(6), 

8 

we  find  that  x  »  1.  for  s  plane  wave,  l.e.  <t  >  ■  1. 

p  I*  •  P*T 

(b)  Bean  waves  and  spherical  waves 

The  atmospheric  transmittance  for  a  beam  wave  Is  affected  by  the 
aperture  size  of  the  transmitter  and  receiver  as  well  as  the  coherence 
of  the  source.  Incorporating  these  parameters  and  using  the  extended 

O 

Huygens-Fresnel  principle,  Vang  and  Plonus9  derived  the  intensity  and 
correlation  function  of  the  received  field  for  a  partially  coherent 
beam  wave  propagating  through  turbulence.  The  intensity  transmittance 
due  to  turbulence  and  source  incoherence  can  then  be  obtained  from 
these  derivations,  l.e. 


<t  >  .  ..i  ♦  5?  +  sL. 

Tr  ill 

T  1  +  C  +  f  +  * 


(7) 


2 

where  C  •  a  /p  ,  C  "  2a  /p  ,  f  ■  ka  /L  Is  the  Fresnel  number  of  the 

s  s  so  s 

source,  a^  Is  the  source  size,  p^  is  source  coherence  length,  Pq  is 
the  coherence  length  of  turbulence,  L  is  the  transmit ter-recelver 
distance  and  k  ■  2nA.  If  a  finite  aperture  receiver  of  radius  R  is 
used,  following  the  definition  of  power  transmittance  <t^  T>,  Eq.(6), 
we  have 


1  -  exp[-(R/a  c)  1 

<T  T>  - . - . Y 

P,T  1  -  exp(-lR/asCa  -  0)}Z] 


(8) 


where 


C2  -  (I  -  l)2  +  (1  +  52  +  C2)/f2 


(9) 


The  detailed  derivation  of  the  intensity  and  power  transmittance  for 
the  various  limiting  cases  has  been  given  in  earlier  reports. 

Here,  we  show  some  brief  results. 

(1)  For  a  collimated  beam  (F  ♦  •)  and  small  receiver  size 
(R  ♦  0)  we  have 


<t  >  h  <v 

P.T  c2  T 


(10) 


The  power  transmittance  is  the  same  as  the  intensity  transmittance 
because  the  received  fields  in  the  small  receiver  area  are  affected  by 
turbulence  uniformly. 

(ii)  Spherical  wave  (a^  >  o) 


<V  -  <Vr>  ■  1 


(11) 


Like  in  the  case  of  the  plane  wave,  a  spherical  wave  is  affected  by 
turbulence  uniformly  throughout  the  entire  receiver  area.  Since  no 
energy  .is  lost  due  to  turbulence  scattering,  turbulence  does  not 
attenuate  the  atmospheric  transmittance  (Intensity  or  power)  for 
spherical  waves. 

(ill)  Incoherent  sources  (p  ♦  --) 

S  Zlt 

<  T  >  -  <T  >  “  1  (12) 

T  P,T 

An  Incoherent  beam  wave  source,  acts  like  a  spherical  wave  source;  it 
radiates  waves  in  all  cirectioi >»,  though  its  field  is  not  coherent. 
Hence,  turbulence  scatters  '  wave  uniformly  and  does  not  give  rise 


t 


> 


Co  Any  Attenuation  in  the  atmospheric  transmittance  for  an  incoherent 
source.  Also,  we  have  shown  that  a  completely  Incoherent  source  (p 

s 

■  0)  does  not  radiate^. 

(iv)  Coherent  source  (p  <■) 

s 


x  1  +  f 

<T_>  - - - - - 

T  1  +  f2  +  52 


(13) 


It  is  interesting  to  note  that  a  completely  coherent  wave  is  subject 

to  the  most  serious  attenuation  due  to  turbulence. 

Eq.(7)  expresses  the  atmospheric  transmittance  due  to  turbulence 

for  a  partially  coherent  source  in  the  turbulent  atmosphere.  When  the 

field  coherence  length  p  is  larger  than  either  the  source  aperture 

o 

alse  a  or  the  Fresnel  tone  /  XL,  the  transmittance  due  to  turbulence 

s 

ia  approximately  unity.  This  is  usually  true  for  the  weak  turbulence 
case.  The  behavior  of  the  power  transmittance  <Tp  T>  is  quite  similar 
to  the  intensity  transmittance  <r^>  when  a  small  aperture  receiver  is 
used.  If  the  receiver  else  is  large  enough  to  collect  all  the 
scattered  field  due  to  turbulence,  i.e.  I  >>  a^c,  no  attenuation  of 
power  occurs.  Therefore,  the  power  transmittance  approaches  the 
constant  unity  »»*♦•.  (See  !q.(8)) 

From  the  above  discussion,  we  conclude  that  the  atmospheric 
transmittance  due  to  turbulence  la  definitely  unity  for  plane  waves  as 
well  as  for  spherical  waves  and  approaches  the  value  of  unity  for  most 
cases  of  beam  waves  except  for  a  coherent  beam  wave  source  in  strong 
turbulencelO .  Hence,  we  have  decided  not  to  incorporate  an 
attenuation  ictor  due  to  turbulence  into  the  calculation  of 


atmospheri 


nsmittance  in  Lowtran. 


III.  Transml ttance  fluctuation  due  to  turbulence 

The  most  serious  effect  of  turbulence  on  wave  propagation  in  the 

turbulent  atmosphere  is  the  fluctuation  of  the  received  field.  We 

2 

defined  the  scintillation  index  m  as  the  normalized  intensity 
2 

variance  , 


c J2  <(I  -  <I>)2> 

<I>2  <I>2 


<I2> 

<I>2 


The  power  scintillation  index  is  given  by 


a2  <(P  -  <P>)2> 

a  2  -  -B-  - . - 

p  <P>2  <P>2 

where 

P  -  /j.  I<£)  d2£ 


<PZ> 


<P>2 


(14) 


(15) 


(16) 


is  the  wave  power  received  by  the  receiver  aperture  E.  We  can  then 

relate  the  scintillation  index  to  the  deviation  of  the  atmospheric 

transmittance  fluctuation.  The  detailed  derivations  of  intensity  and 

power  transmittance  bounds  are  shown  in  Appendix  4A  and  4B  in  the 

Fourth  report  .  From  this  report  we  obtain  that 

f..  ■  <r>  (1  ±  m)  «  t.  (I  ±  m)  (17) 

l  L 

where  xu  and  are  the  upper  and  lower  bounds  of  the  transmittance, 
respectively.  And, 

TP  *  <TP>  (1  1  mP)  (l8) 

u 
l 

Mote  that  the  upper  and  lower  bounds  of  transmittance  are  not  the 


exact  bounds  as  some  measured  transmittance  may  h*  out  of  the  bounds. 


However,  for  many  samples,  we  expect  that  most  measured  transml ttances 
will  fall  inside  these  bounds.  Since  we  concluded,  in  section  II, 


that  the  transmittance  attenuation  due  to  turbulence  scattering  is  not 
significant,  the  only  involvement  of  turbulence  in  Eqs.  (17)  and  (18) 


is  in  the  scintillation  index  m2  and  m2.  The  new  subroutines  VRANI 

P 


and  SINTL  are  to  calculate  m2  and  m2  for  plane  waves  and  beam  waves. 

P  ’ 


respectively. 

(a)  Plane  waves 

(i)  Point  receiver 

Consider  a  plane  wave  U  propagating  through  the  turbulent 
medium  represented  by 

0  -  e  x  +  1  S  (19) 

where  x  and  S  present  the  random  log-amplitude  and  random  phase  due  to 
turbulence,  respectively.  Assuming  a  Gaussian  probability 
distribution  for  x  *  the  averaged  intensity  and  variance  of  Intensity 
can  then  be  stated  as 


2X 


<I>  -  <0  •  U  >  -  <e  >  -  e 


2<X>  +  2o 


(20) 


<1 


2>  -  <0  »  U  •  U  •  U*> 


4  <X>  +  go 


(21) 


Substituting  Eqs. (20)  and  (21)  into  Eq.(14),  the  scintillation  index 
for  plane  waves  is 


(22) 


2 

The  variance  of  log-amplitude,  o 


has  be^n  found  by  Rytov’s  method 


(> 


10 


However,  It  is  only  valid  for  weak  turbulence  when  applied  in  Eq.(22). 

Experimental  data  indicates  that  m2  (i.e.  a2  )  saturates  toward  the 
6  N 

value  of  unity  .  In  recent  years,  theoretical  work  to  prove  chat  the 

variance  of  intensity  saturates  to  a  constant  of  unity  was 

performed*®.  Avoiding  complex  mathematics  and  hoping  to  get  a  model 

which  is  sufficiently  accurate  under  weak  and  strong  turbulence 

conditions,  we  relate  the  variance  of  intensity  and  log-amplitude  by 
-2a2 

m  -  1  -  e  x  (23) 


For  small  values  of  a^,  m  ■  2a^  which  agrees  with  Eq.(22).  For  large 

m  »  1,  which  agrees  with  the  saturation  condition.  Using  Ref. 

(7),  the  variance  of  log-amplitude  as  found  by  Rytov's  method  is  given 


a2  -  0.563  k7/6  /  L  C2(tj)  <L  -  U>5/6  dU  (2< 

X  on 

where  k  *■  2it/\  is  the  wavenumber. 

C2  is  the  structure  constant  of  turbulence, 
n 

For  a  homogeneous  medium,  C2  is  constant  along  the  path,  and  Eq.(24) 

n 

can  be  rewritten  as 


2  n  „  A  ,7'6  tU/6 

a  -  0.31  C  k  L 

X  n 


2  .  11 
A  model  of  for  the  earth  s  atmosphere  is  given  by  Hufnagel,  et 

al.  We  will  show  and  modify  this  model  to  fit  in  Lowtran  later. 

(ii)  Finite  aperture  receiver 

For  a  finite  size  of  receiver  with  radius  R  ,  the  average 


received  power  for  Che  plane  wave  case  can  then  be  obtained  from 


Eq.  ( 16) , 

<P>  -  Jj.  <I>  d2  £  -  <X>  HR2  (26) 

and  similarly  Che  mean-square  received  power  Is 

<P2>  -  /j/j.  <H>  d2  £j  d2  £2  •  (27) 

Assuming  chat  x  is  Gaussian,  we  substitute  Eqs.(19),  (26)  and  (27) 
into  Eq.(15)  and  obtain  m2  , 


•  “2  /”  -  »  [eo.-'tlg)  -  (fj)  / 1  -  ef-]  pap  (28) 


HR  o 


4R 


where 


B^(p)  S  <x(£>X<£  +  £)> 


(29) 


The  averaging  factor  G(R)  Is  defined  as 


2,6 


G(R)  ■  m2/m2 


(30) 


From  Eqs.  (22)  and  (28),  we  obtain  the  following  expression  for  G(R) , 


G(R)  -  ---  J2R  b^p)  [co,‘l(|-)  -  e.  /T7Z]  pdp 


hR2  o 


(31) 


4R^ 


where 


bx(p  ) 


(32) 


4  c 


e  X  -  l 


Using  Rytov's  Method,  the  Kolmogorov  spectrum  and  a  locally 
homogeneous  medium,  the  correlation  function  of  log-amplitude  can  be 


found  from  Ref.  (6)  under  the  condition  L  X  »  1q, 


B  X(P)  -  0.033n2(-  T  N-  jhl  k2  JL  C2^)  dr,  J  (p,h)  (33 

o  on 

where 

j(P.n)  -  +  ---£:-3-)5/6lfl(-  \.  i. . e- . ) 

"■  «(.}  +  .i<t  :.!>1 

K2 

m 

-(~)5/6,  F,  f-  f,  1,  *)  (34 

K2  6 

m 


X  -  K2p2/4 

m 

K  "  5.92/fo  {35 

is  the  inner  scale  of  turbulence 
jF^(a,b,x)  is  the  degenerate  hypergeometric  function. 

To  calculate  the  power  transmittance  deviation  in  subroutine  VRANI,  w( 
first  obtain  the  G(R)  factor  and  use  the  modified  intensity  deviation 
Eq.(23) ,  that  is 


mp  -  m  •  G(R)  (36! 

(b)  Beam  waves  and  spherical  waves 

The  extended  Huygens-Fresnel  principle  can  be  used  to  obtain  the 
receiver  field  for  beam  wave  or  spherical  wave  sources  in  the 
turbulent  medium: 


where  x(s,p)  *  X  +  iS  is  the, random  pertubetion  due  to  turbulence. 
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From  Eq.(37),  wc  can  derive  the  expressions  of  the  received  Intensity 


and  intensity-correlation  functions  that  are  needed  to  obtain  the 
Intensity  and  power  scintillation  index.  For  a  practical  system,  the 
source  is  not  necessarily  coherent.  Thus,  there  could  exist  a  random 
part  of  u(s)  in  Eq.(37).  However,  a  slow-response  time  (narrow 
bandwidth)  receiver  can  smooth  out  some  fluctuations  due  to  source 
randomness.  Hence,  the  intensity  and  power  scintillation  index  are 
affected  ,  in  addition  to  turbulence,  by  source  incoherence  and 
receiver  response  tine  (bandwidth). 

The  intensity-correlation  function,  which  includes  the  effects  of 
source  incoherence  and  random  medium  due  to  turbulence,  has  been 
derived  for  a  partially  coherent  beam  wave  source  in  turbulent 
medium* ,  as 


®  -  <I(£1)I(£2)> 

"  /••'/  d2V2l2<,2i3d2V*<Il»l2»f3»*A>,r4<£l*i2»i3'#45£l»£2> 

•  «*p{-  *  £2)  +  £2  *  <*3  "  *4^1  +  ZC^*1  +  *2  +  *3  +  *4^ 

(38) 

where 

y4(il'l2»l3*iA)  "  <us(il)us<i2)tts(l3)us(i4)>s  (39> 

is  the  fourth-order  source  coherence  function  and 


m 


is  the  fourth-order  spherical  wave  coherence  function  in  the  turbulent 
medium.  The  bracket  subscripts  s  and  m  denote  the  ensemble  averages 
over  the  statistics  of  source  and  turbulent  medium,  respectively. 

(i)  Source  statistics  and  response  time  of  receiver 
For  a  receiver  with  response  time  larger  than  source  coherent 
time,  F*  in  Eq.(39)  reduces  to  a  product  of  second-order  spherical 
wave  correlation  functions^ ,  13  i,e. 

K  "  <us(±l)us  (i2)><iU(ll)tta(i2)>  (41> 

Eq.(41)  makes  the  mathematics  much  simpler  and  correctly  gives  zero 

scintillation  in  vacuum.  On  the  other  hand  if  the  response  time  is 

smaller,  the  fluctuations  due  to  the  source  are  not  smoothed  out.  The 

following  derivation  is  to  model  mathematically  the  source  coherence 

properties  and  obtain  a  suitable  expansion  of  Eq.(39). 

For  a  partially  (spatially)  coherent  source,  u  (a)  can  be 

s 

expressed  by  the  product  of  the  deterministic  radiation  distribution 

factor  u  ,(s)  and  the  random  coherence  factor  u  (s).  Let  u  ,(s)  be  a 
so  —  sr  —  sa  — 

source  distribution  such  as  that  of  a  fundamental-mode  laser: 

usd(»)  ■  As  expf  -  f --s  +  — )s2] 

“  2  or  2F 

s 

and  let  the  random  part  of  the  source  field  be 

usr(s)  -  e1*^  (A3'i 

We  model  the  random  phase  of  the  source  field  as^ 

s)  ■  a  +  b  •  s 


UO 
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1, 

( 


where  e  end  b  ere  e  random  shift  end  e  random  tilt  vector  of  the 
random  phase.  Assuming  the  distributions  of  a  and  b  are  Gaussian  with 
zero  mean,  we  have  two  kinds  of  statistics,  namely  Gaussian  phase 
statistics  and  pair-correlated  field  statistics  to  apply  to  F* . 

In  Ref.  (4),  we  have  shown  that  pair-correlated  statistics 
give  the  better  results.  The  source  coherence  function  obtained 
bv  these  statistics  can  be  stated  as. 


rs,  .  _  .4  r  1  f  2,  2  2>  jk,  2  2^2  2.-, 

r4(s1*i j s'i**  As  exPi  — 2»*i  +  *2  *j  +  s4)  *  2F^*1  ’  s2  +  *3  "  3A  - 


(exp{---  fUj  -  £2)2  +  <s3  -  *^>2] } 


+  «xp{-i-F(£l  -  nj2  +  (s3  -  «^)2]} 


fa 


+  exp {-  4<^  -  --;[(,*!  +  S3)2  +  («2  +  s^)2  ]} 
4ps 

-  2exp{  -  2oJ  -  ---(s*  +  »2  +  *3  + 


where 

p2  -l/<b2>  (46) 

o2  is  the  variance  of  the  random  shift  a,  a  is  the  beam  radius  and  F 

a  • 

2 

is  the  radius  of  curvature  of  beam  wavefront,  a  and  p  are  measures 

A  S 

of  the  degree  of  coherence.  As  ->  •  .  ->  0,  we  consider  the 

source  coherent;  if  p  ->  0  or/and  a2  ->  ®  ,  the  source  is  incoherent. 


(ii)  Atmospheric  turbulence 


In  weakly  turbulent  media,  we  can  assume  that  the  random 
pertubation  ?  is  Gaussian,  i.e.  the  log-normal  field  assumption  is 
valid.  The  fourth-order  spherical  wave  coherence  function  can  then 
be  expressed  by  the  structure  functions  and  correlation  functions  as'* 

h  “  exp^  2  *±1  '  —2 *  0>  ‘  2  D(^l  ‘  —4 *  Pd>  _  2  D(^2  '  -3’  V 

\  D<-3  *  0)  +  f  d(±2  '  14*  pd*  +  f  D(Il  ‘  —3 *  pdJ 

+  2BX(i2  *  14*  £d)  +  2Bx(±l  "  i3’  Pd> 

+  1  DxS^£2  "  14*  Pd^  “  1  DxS^il  "  ±3*  pd ^  1 

(47) 

where  ■  pf  -  p^  «nd  D,  ,  D^g  are  the  wave  structure  function, 
log-amplitude  correlation  function  and  log-amplitude  phase  structure 

15 

function,  respectively.  The  two-wave  structure  functions  are  known. 
1/2 

Hence,  for  (\L)  »  I  s .  |  »  1.  ,and  by  use  of  the  quadratic 

<1  0 

approximation,  we  can  find  the  wave  structure  functions, 


Hid*pd) "  i(sd  +  id  *  J*d  +  pd) 

po 

^s<id*£d)  -:i'(sd  +  »d  -fd  +  p2d> 


xs 


where 


--  -  1.575  kl2/5  L  r/!rdTi(L  -  ti)5/3  <£  (n)]  6/5 

-l-  -  0.234  k13'6  l’1^6  /Ldr,— Ss!<!!>!!? 


xs 


O'  •  .  ,,1/6 

[t|(  L  -  h)l 


(48) 

(49) 

(50) 

(51) 


The  use  of  the  quadratic  approximation  for  the  structure  functions 
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does  not  imply  that  we  arc  limited  to  the  case  of  tilt-only 

medium1^’*7  because  in  the  expansion  of  F^  terms  other  than  phase-tilt 

14 

terms  are  present  and  which  are  retained.  Fante  has  Introduced  a 
useful  log-amplitude  correlation  function  as 


VvV  ■  v  • 


*  "5(ad  +  ij  •  Sd  +  Pd) 


where 

a2  -  0.225  b7/6/JdT,  Cn2(Ti)(L  -  r\)5/ 6  (! 

*s 

is  the  variance  of  log-amplitude  for  spherical  waves*.  To  obtain  a 

2  2 

closed  form  result  for  m  and  m  ,  we  should  further  approximate 

P 

£q.(52)  as 


B  *  Oj  U  -  ^  •  Pj  +  P2)] 


Note  that  only  the  structure  constant  C2(ti)  contained  in  Eqs.  (50), 


(51)  and  (53)  characterizes  turbulence  properties.  Therefore,  once  we 

?  1  1  2 
know  C-  along  the  propagation  path,  we  can  obtain  -3,  -3-  and  a  for 

Po  PXS 

both  homogeneous  (horizontal  path)  and  inhomogeneous  (slant  path) 
turbulent  media. 

2 

The  structure  constant  C  has  been  measured  and  modeled  for  the 

n 

11 

earth's  atmosphere  by  Hufnagel,  et.  al.  We  modify  it  to  fit  Lowtran 


C>i  - 
n 


4.2  X  10'Uh-2/3  exp(-h/320)  (h  >  10  m) 


8.77  X  10 


(h  <  10  m) 
(h  >s  100  Km) 


0 
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where  h  is  the  altitude  in  meters. 

(iii)  Intensity  and  power  scintillation 

A  step-function  receiver  makes  the  mathematics  complicated  such 
that  a  closed-form  of  power  scintillation  cannot  be  obtained. 
Therefore,  we  integrate  the  given  intensity  and  intensity-correlation 

a 

function  over  the  receiver  aperture  weighted  by  a  Gaussian  function, 
which  allows  us  to  relax  the  integration  limits  to  infinity  and  obtain 
closed-form  results.  Using  Eqs.  (14),  (15),  (37),  (38),  (45)  and 
(47),  the  intensity  and  power  scintillation  index  for  a  partially 
coherent  beam  wave  source  in  turbulence  have  been  obtained.2  That  is, 


f  -1-  +  +  --  +  (— -)2]  exp{f  — ]p2 } 

i4  '4a2  4p2  p2  2  [a2  A2  +  1-  +  ‘i-  +  I-J 


a2  u2  2«2  ,2  2 


«xp{-  5-f- ] 
UL* 


Ixyz- 


a  *  2 

„2 ,  +  i.  +  ( -|-)2]  r-1-  + 


ka2  2 


4a2  4p2 

s  s 


a2fi  +  (  — )2  +  {-— )2  +  fa2A)2; 


m 


••  ,s 

V 


in'  .V 
k‘.V 


v  • 


i 

K' 


>y 


>>< 

A 

■v- 


>>■> 


i 

* 


.'V 

I 


i 


i 


i 


K> 


r.-‘ 


0  ■  D  ♦  -- 


4p‘ 


(A  -  -J-)‘ 

P_e 


2  “  o  +  "  ifi — 


(58) 


Instead  of 


ng  B 


'  °X  “  j  ®X*  "*  *“v*  uaed  ^9*  for  the 

log-amplitude  correlation  function.  Therefore,  we  substitute  ( --  - 


2d* 

1  *g 

-j)  by  (----)  for  all  terms  in  Eq.  (59)  which  can  be  obtained  from 

P  P2 

X  o 


Ref.  (2).  This  substitution  gives  the  better  results,  even  though  the 

.16 


application  range  for  the  approximation  is  still  limited 

By  letting  a  ->  0,  we  can  obtain  the  spherical  wave  results.  The 
intensity  and  power  transmittance  deviations  for  beam  waves  and 

spherical  waves  can  then  be  calculated  by  the  SINTL  computer  code 
using  Eqs.  (56)  and  (57). 


.  J  V-  - 


•>? 


IV.  Subroutine  VRANI 

The  subroutine  VRANI  is  designed  to  calculate,  by  using  Eqs .  (23) 

-  (36),  the  plane  wave  intensity  and  power  variance. 

For  the  program,  we  changed  the  integration  in  all  formulas  to 

the  summation  form.  From  Eq.(55),  we  find  that  varies  rapidly  when 

height  h  is  small.  By  trading  off  calculation  time  and  precision  of 
2 

C  ,  we  choose  the  subintervals  in  the  summation  as, 
n  7 


4hlJ  ■ 

20  m 

bu 

< 

25  km 

“ij  ■ 

100  m 

25 

km  < 

hu 

< 

50  km 

“«1- 

400  m 

50 

km  < 

-u 

< 

100  km 

Eq.(24)  can  then  be  rewritten  as 


£L 

0.56  Z  Z  C2(h  )<L  .  Ll.)5/6  C-— Ah 

i  j  n  1J  -  hj.j  ij 


where  h^  is  the  height  which  corresponds  to  the  calculated  points  on 
the  path,  "i"  is  the  layer  index,  "J"  is  the  sub  index  of  each  layer, 


L  is  the  total  path  length,  is  the  path  range  from  the  transmitter 
to  the  point  calculated  and  is  the  path  range  for  each  layer 
passed. 

Since  the  structure  constant  depends  on  h,  the  calculation  for 
horizontal  and  slant  paths  are  different. 

(a)  Horizontal  path 

2 

The  intensity  variance  m  can  be  obtained  by  Eqs.  (23)  and  (25) 

after  the  constant  is  specified.  B^(p)  can  ^e  found  with  the 

condition  that  <<  \  I.  by  using  Rytov's  Method  and  Kolmogorov 
o' 

H  1  7 

spec  trum-’*  , 


1 


“Tx)S7*V75 

bx(p)  -  1  -  2.36f-£")5/6  + 

-  0.0242  ("-)’7/6 

where  X  is  the  inner  scale  of  turbulence  and  is  assumed  to  be  3  ram  in 
o 

our  subroutine.  The  power  variance  can  then  be  obtained  by  Eqs.  (31) 
and  (36). 

(b)  Slant  path 

Eq.  (60)  can  be  used  to  calculate  the  intensity  variance  for  the 
slant  path  cases.  To  obtain  results  for  power  variance,  Eqs.  (33)  and 
(34)  have  to  be  approximated  so  that  the  integration  can  be  changed  to 
a  summation  from.  The  calculation  for  power  variance  can  then  be 
performed.  The  approximations  are 

1  -  0.8333x  -,0347x2  -  0.0045X3  |x|<  6.5 

iPj-  1,  *)  -  (63) 

1.0627  (  -  x)5/6  | x | >  6.5  and  (Re  x  <  0) 


P  V  V  A. 


-  0.024(-£-)  <  <  X  <  </TT  (62) 


/  X.L  <  <  p~ 


and 

0  (y  >  6.5,  *  >  6.5) 

(I:*T)5/6(o.259  +  0.805y  +  0.009y2  -  0.0043y31 


l.C^Sy3'6 


(y  <  6.5,  x  >  6.5) 


J(0,P)  - 


(64) 


(t:D)5/6(0>259  +  o.805y  +  0.009y2  -  0.0043y3) 


(L)5/6(1  +  o.8333x 

Iv 

m 


-  0.034 7x2  +  0.0045x3) 


(y  <  6.5,  x  <  6.5) 


where 


x  »  -  K2  p2/4 
m 


P2k 

y " 


(65) 

(66) 


The  approximations  have  been  checked  and  the  total  error  is  under  10Z. 

L-h  l 

Furthermore,  we  also  assume  that  -r-  >  >  -r  which  means  that  we 

KZ 


m 

neglect  a  very  small  part  of  turbulence  near  the  receiver. 

The  subroutine  VRANI  is  called  by  another  subroutine  TRANS  in 
Lowtran,  which  calculates  the  atmospheric  transmittance  due  to 
absorption  and  scattering  of  molecules  and  aerosal.  The  program  is 
listed  in  Appendix  B  and  the  definitions  of  symbols  and  variables  are 
shown  in  Appendix  A. 

For  comparision,  we  used  the  1962  U.S.  standard  atmospheric  model 
the  rural  aerosal  model  with  visual  range  of  23km  and  a  5km 
propagation  distance  to  plot  all  the  figures.  For  horizontal  paths  an 
altitude  of  400m  is  used  for  all  figures.  For  an  upward  path,  the 
altitudes  of  transmitters  and  receivers  are  set  to  800m  and  400m, 
repectively.  For  a  downward  path,  the  altitudes  of  transmitter  and 


receiver  are  reversed. 


Figs.  1-3  show  the  plane  wave  intensity  transmittance  for 
horizontal,  downward  and  upward  paths,  respectively.  Figs.  4-6  give 
the  power  transmittance  for  the  seme  condition  except  that  a  10cm 
receiver  is  used.  The  latter  figures  display  a  smaller  transmittance 
deviation  (scintillation).  This  is  due  to  the  receiver  aperture 
averaging  effect. 


V. 


Subroutine  SINTL 


The  subroutine  SINTL  calculates  the  intensity  and  power  variance 
for  partially  coherent  beam  wave  and  spherical  wave  sources. 


(a)  Horizontal  path 
2 

C  is  constant  along  the  path  for  the  horizontal  case.  Cqs. 
n 

(50),  (51)  and  (53)  can  then  be  rewritten  as 


-i-  -  (0.546  C2  k2  L)6/5 
-2  n 


(67) 


0.114C2 

n 

-  0.124 


k13/6  l5/6 
c2  *7/6  l11/6 


(68) 

(69) 


2  11 
After  the  constant  is  calculated  from  Eq.  (55),  -j-,  -2”  and 

2  n  po  pxs 

can  be  obtained  immediately  from  the  above  equations.  Eq.  (56) 

s 

will  then  give  the  result  of  the  Intensity  scintillation  for  a 
partially  coherent  beam  wave  source. 

(b)  Slant  path 

Changing  the  Integrations  in  Eqs.  (50),  (51)  and  (53)  into 


summation  form,  like  that  of  Eq.  (60),  we  can  obtain  the  results  for 


1  1  .2 

,  -n -  ana  c 

7  ^  v/ 


The  power  scintillation  (variance)  can  then  be 


o  XS 

calculated  by  the  closed-form  formula  of  Eq.  (57). 

SINTL  is  also  called  by  the  subroutine  TRANS  when  the  transmitter 


size  Is  finite.  The  program  list  is  shown  in  Appendix  C. 

figs.  7-9  show  the  intensity  transmittance  for  partially 


coherent  been  wave  sources  in  horizontal,  downward  and  upward  paths, 
respectively.  For  the  finite  receiver  with  radius  of  10cm,  Figs.  10 
12  give  the  power  transmittance.  Again,  the  receiver  aperture 
averaging  effect  can  be  found  by  comparing  the  intensity  and  power 
transmittance  deviations.  The  results  of  spherical  wave  cases  are 
shown  in  Figures  13-15  and  Figures  16  -  18  for  point  receivers  and 
finite  aperture  recievers,  respectively. 


APPENDIX  A 


ANGLE 

BI 

BL 

BX 

CN2 

DD 

DH 

DO 

DS 

DSW 

DT 

DZW 

FR 

GAA 

GD 

HMIN 

HW 

HI 

H2 

IV 


Symbols  and  Defintions  for  Subroutines  VRANI  and  SINTL 

Initial  zenith  angle  in  degree 
Covariance  of  Intensity 

Log-amplitude  covariance  normalized  by  variance 

Covariance  of  log-amplitude 

C^  -  structure  constant  of  turbulence 
n 

the  ratio  of  the  distance  from  point  calculated  to 
receiver  over  total  path  length. 

Height  interval  of  slant  path  integration 
Distance  from  point  calculated  to  transmitter 
Distance  from  point  calculated  to  receiver 
Path  length  in  a  layer 

Same  as  DS,  especially  used  in  the  downward  long  path 

calculation 

Height  for  a  layer 

Fresnel  zone  in  neter  (m) 

,  variance  of  random  shift  for  a  partially  coherent  source 
It  is  assumed  to  be  zero  for  a  temporal  coherent  source. 
Aperture  averaging  factor 
The  minimum  height  of  a  downward  path 
Height  corresponding  to  the  point  calculated 
Height  of  receiver  (and  transmitter  for  horizontal  path) 
Height  of  transmitter 
Wavenumber  in  cm  * 


JMIN  The  layer  index  of  the  minimum  height  for  a  downward  path 


PVR 


Power  variance 


APPENDIX  B 


31530 

31540  C 

31550  C 

31560  C 

31570  C 

31580  C 

31590 

31600 

31610 

31620 

31630 

31640 

31650 

31660 

31670 

31680 

31690 

31700 

31710 

31720 

31730 

31740 

31750 

31760 

31770 

31780 

31790 

31800  C 

31810  C 

31820  C 

31830 

31840 

31850 

31860 

31870 

31880 

31890 

31900 

31910 

31920 

31930 

31940 

31950 

31960 


SUBROUTINE  VRANI(IV) 

THIS  SUBROUTINE  IS  TO  CALCULATE  THE  VARIANCE  OF  INTENSITY 
DUE  TO  TURBULENCE  AND  THE  CALCULATED  STANDARD  DEVIATION  CAN  BE 
USED  TO  DEFINE  HIGH  BOUND  AND  LOW  BOUND  OF  TRANSMITTANCE 

COMMON  /CARD1/  MODEL, IHAZE , ITYPE ,LEN, JP ,IM,Ml ,M2 ,M3 ,ML .IEMISS ,R0 
1  ,TB0UND .ISEASN ,IVULCN , VIS 

COMMON  /CARD2/  HI ,H2, ANGLE, RANGE, BETA, HMIN, RE, TT,RR 
COMMON  /CARD3/  VI ,V2 ,DV ,AVM,CO, CW,W( 15) ,E (15) ,CA,PI 
COMMON  /CNTRL/  LENST.KMAX , M, IJ , J1 , J2 , JMIN .JEXTRA, IL,IKMAX,NLL,NP1 
1 ,IFIND ,NL , IKLO 

COMMON  /WANG/  K2 ,DSW(34 ) ,DZW(34 ) ,XWl (34  )  ,XW2  (34) 

COMMON  /VRAN/  VRI 
DIMENSION  WS(34) 

HW-H1* 1000.0 

WRANGE-RANGE*1000.0 

WH2-H2*1000.0 

VR-0.0 

CN2-0.0 

P VR-0.0 

WL-0.01/IV 

FR-(WL*WRANGE )**0. 5 
WK-IV*100 . *2 . *P I 
VK-WK**(7./6. ) 

W0-9.E-6/5. 910**2. 

IF (ITYPE. NE.l)  GO  TO  20 

VARIANCE  CALCULATION  FOR  HORIZONTAL  PATH 

CN2-4 . 2E-14*HW** (-2 . /3 . )*EXP(-HW/320. 0 ) 

IF(H1.GT. 100.0)  CN2-0.0 
IF(HW.LE.IO.O)  CN2-8.77E-15 
VR-0.31*CN2*WRANGE**(ll./6. ) 

VR-VR*VK 

VRI- 1 . -EXP ( -2 . *VR*  *0 . 5 ) 

IF(RR.LT. 0.001)  GO  TO  91 
DO  18  1-1,100 
Y-0.01*I 
DY-RR*2.*Y 

IF (DY.GE. 0.003)  GO  TO  11 

BL-1 .-12 .3*DY**2 .0/ (FR**(5 . /3 . )*0 .003**(1 . / 3 . ) ) 

GO  TO  17 

11  XI-WK*DY**2 .O/WRANGE 


31970  IF (DY.GE.FR)  GO  TO  12 

31980  BL-l.-2.36*XI**(5./6.)+1.71*XI-0.024*XI**2.0 

31990  GO  TO  17 

32000  12  BL— 0.0242* (XI/4 . )**(-7. /6. ) 

32010  17  CONTINUE 

32020  BX«BL*VR 

32030  BI-EXP(4.0*BX)-1. 

32040  P VR-PVR+BI * ( ACOS ( Y )-Y* ( 1 . -Y  **2 . ) **0 . 5  )  *Y  *0 . 0 1 

320S0  18  CONTINUE 

32060  PVR-PVR*16./Pl 

32070  WV-EXP(4.0*VR)-1. 

32080  GD-PVR/WV 

32090  VRI-VRI*GD**0 .5 

32100  GO  TO  91 

32110  20  IF (ANGLE. GT. 90.0)  GO  TO  37 

32120  C 

32130  C  VARIANCE  CALCULATION  FOR  UPWARD  PATH 

32140  C 

32150  DO  35  K-1,50 

32160  DS-0.1 

32170  BX-0.0 

32180  Y-0.02*K 

32190  DY-RR*2 . *Y 

32200  RF-DY**2./4. 

32210  VR-0.0 

32220  Q-RF/WO 

32230  JK>  34  I-J1.J2 

32240  WS(I)-DSW(I)/DZW(I) 

32250  IF(XW1(I).LE.25.0)  GO  TO  21 

32260  IF(XWl(I) .LE.50.0)  GO  TO  22 

32270  DH-400.0 

32280  GO  TO  23 

32290  21  DH-20.0 

32300  GO  TO  23 

32310  22  DH-100.0 

32320  23  WIK-(XW2(I)-XW1(I))*1000.0/DH 

32330  IK-WIK 

32340  HW-XW1(I)*1000.0 

32350  DO  33  J-l.IK 

•  32360  HW-HW+DH 

32370  IF(HW.GE. 99600.0. AND. ITYPE.EQ. 3)  XW2 (I )- 100000.0 

32380  IF (HW.GE. 99600.0)  GO  TO  34 

32390  CN2-4 . 2E-14*HW**(-2 . /3 . )*EXP(-HW/320.0 ) 

32400  IF(HW.GT. 100000. )  CH2-0.0 


32410  IF (HW.LE. 10 . )  CN2-8.77E-15 

32420  VR-VR+0.56*CN2*DS**(5./6.)*DH*WS(I) 

32430  DD-DS/WK 

32440  R-RF/DD 

32450  IF(DD.LT.WO)  GO  TO  27 

32460  IF(R.GE.6.5.AND.Q.GE.6.5)  GO  TO  27 

32470  IF(R.LT.6.5.AND.Q.LT.6.5)  GO  TO  28 

32480  G«DD**(5./6. )*(0. 259+0. 805*R+0.009*R**2.-0.0043*R**3. )-l .063*RF**( 

32490  *5./6.  ) 

32500  IF(G.LE.O.O)  G-0. 

32510  GO  TO  31 

32520  27  G-0.0 

32530  GO  TO  31 

32540  28  G-DD**(5./6. )*(0. 259+0. 805*R+0.009*R**2.-0.0043*R**3. )-W0**(5./6. ) 

32550  **(1.+0.8333*Q-0.0347*Q**2.+0.0045*Q**3. ) 

32560  IF(G.LE.O.O)  G-0.0 

32570  31  BX«BX+G*CN2*DH*WS(I) 

32580  DS-DS+DH*WS ( I ) 

32590  33  CONTINUE 

32600  DQ-WRANGE-DS 

32610  IF (DQ.LE .0 . )  DS-DS~DH*WS (I ) 

32620  VR-VR+0 . 56*CN2  *DS** (5 . /6 . ) * (XW2 ( I )*1000 . -HW ) *WS (I ) 

32630  BX-BX+CN2 *G* (XW2 ( I )*1000 . -HW )*WS ( I ) 

32640  DS-DS+(XW2 (I)*1000.0-HW)*WS (I ) 

32650  34  CONTINUE 

32660  IF(RR.LT. 0.001)  GO  TO  36 

32670  BX«BX*2.117*WK**2. 

32680  BI-EXP(4.*BX)-1. 

32690  PVR«PVR+BI*(ACOS(Y)-Y*(1.-Y**2.)**0.5)*Y*0.01 

32700  35  CONTINUE 

32710  PVR«PVR*16./PI 

32720  36  VR-VR*VK 

32730  VRI«1.-EXP(-2.*VR**0.5) 

32740  IF(RR.LT. 0.001)  GO  TO  91 

32750  WV-EXP(4.*VR)-l. 

32760  GD-PVR/WV 

32770  VRI-VRI*GD**0.5 

32780  GO  TO  91 

32790  37  CONTINUE 

32800  C 

32810  C  VARIANCE  CALCULATION  FOR  DOWNWARD  PATH 

32820  C 

32830  DO  62  MW-1,50 

32840  DS-0.1 


32850 

32860 

32870 

32880 

32890 

32900 

32910 

32920 

32930 

32940 

32950 

32960 

32970 

32980 

32990 

33000 

33010 

33020 

33030 

33040 

33050 

33060 

33070 

33080 

33090 

33100 

33110 

33120 

33130 

33140 

33150 

33160 

33170 

33180 

33190 

33200 

33210 

33220 

33230 

33240 

33250 

33260 

33270 

33280 


DT-0.1 

BX-0.0 

Y-0.02*MV 

DY-RR*2.*Y 

RF-DY**2./4. 

VR-0.0 

Q-RF/WO 

Ll-Jl 

DO  60  L-l.NL 
WS(L1)-DSW(L1)/DZW(L1) 

IF(XWl(Ll).LE.25.0)  GO  TO  38 
IF(XW1(L1).LE.50.0)  GO  TO  39 
DH-400.0 
GO  TO  40 

38  DH-20.00 
GO  TO  40 

39  DH-100.00 

40  WIK-(XW1(L1  )**XW2 (LI ) )*1000.0/DH 
IK-WIK 

HW-XW1(L1)*1000.0 
DO  57  J-l.IK 
HW-HW-DH 

CJf2-4.2E-14*HW**(-2. /3 .  )*EXP(-HW/320.0) 

1F(HW.GT. 100000.)  CN2-0.G 
IF(HW.LE. 10. )  CN2-8.77E-15 
VR-VR+0.56*CN2*DS**(5 . /6 . )*DH*WS(L1 ) 

DD-DS/WK 

R-RF/DD 

IF (DD.LT.WO)  GO  TO  51 
IF(R.GE.6.5.AND.Q.GE.6.5)  GO  TO  51 
IF(R.LT.6.5.AND.Q.LT.6.5)  GO  TO  52 

Gl“DD**(5./6. )*(0. 259+0 .805*R+0 .009*R**2 .-0.0043*R**3 . 
*(5./6. ) 

IF (G1 .LE.O. )  Gl-O. 

GO  TO  53 

51  Gl-0.0 
GO  TO  53 

52  Gl-DD**(5./6. )*(0. 259+0. 805*R+-0.009*R**2.-0.0043*R**3. 
*)*(1.+0.8333*Q-0.0347*Q**2 .+0.0045*Q**3.) 

IF(Gl.LE.O.O)  Gl-0.0 

53  BX-BX+G1*CN2*DH*WS(L1) 

DS-DS+DH*WS (LI ) 

IF(K2.EQ.O)  GO  TO  57 

IF (K2.EQ.1.AND.WH2.LE.HW)  GO  TO  57 


i-l  .063*RF** 


i-WO**(5./6. 


3-4 


33290  VR»VR+0.56*CN2*(WRANGE-DT)**(5./6. )*DH*WS(Ll) 

33300  DD-(WRANGE-DT)/WRANGE 

33310  R-RF/DD 

33320  IF(DD.LT.WO)  GO  TO  54 

33330  IF(R.GE.6.5.AND.Q.GE.6.5)  GO  TO  54 

33340  IF(R.LT.6.5.AND.Q.LT.6.5)  GO  TO  55 

33350  G2«DD**(5./6. )*(0. 259+0. 805*R+0. 009*R**2 .-0.0043*R**3 . )-l .063*RF** 

33360  *(5./6. ) 

33370  IF (G2 .LE.O . )  G2-0. 

33380  GO  TO  56 

33390  54  G2-0.0 

33400  GO  TO  56 

33410  55  G2-DD**(5 . /6. )* (0.259+0. 805 *R+0. 009 *R**2 ,-0.0043*R**3. )-WO**(5./6. 

33420  *)*(1.+0.8333*Q-0.0347*Q**2.+0.0045*Q**3.) 

33430  IF (G2 .LE.O.O)  G2-0.0 

33440  56  BX-BX+G2*CN2*DH*WS(L1) 

33450  DT-DT+DH*WS(L1) 

33460  57  CONTINUE 

33470  DQ-WRANGE-DS 

33480  IF (DQ. LE.O.)  DS«DS-DH*WS(L1 ) 

33490  VR-VR+O . 56*CN2*DS**( 5 . / 6 . )* (HW-XW2 (LI )*1000 . )*WS (LI ) 

33500  BX*BX+G 1 *CN 2  * (HW-XW2 (Ll )*1000. )*WS (LI ) 

33510  DS-DS+(HW-XW2(L1) *1000.0) 

33520  IF (K2 .EQ .0 )  GO  TO  58 

33530  IF (K2 . EQ . 1 . AND . WH2 . LE . HW )  GO  TO  58 

33540  VR-VR+O . 56*CN2* (WRANGE-DT )** (5 . /6 . )* (HW-XW2 (Ll )*1000 . )*WS (Ll ) 

33550  BX-BX+G2 *CN2 * (HW-XW2 (Ll )*1000. )*WS (Ll  ) 

33560  DT-DT+(HW-XW2(L1)*1000.0)*WS(L1) 

33570  58  CONTINUE 

33580  Ll-Ll-1 

33590  IF (K2 .EQ .0.AND.L1 .LE. J2)  GO  TO  61 

33600  IF (Ll .LE.JMIN.AND.K2 .EQ . 1 )  GO  TO  61 

33610  60  CONTINUE 

33620  61  CONTINUE 

33630  IF(RR.LT. 0.001)  GO  TO  90 

33640  BX-BX*2 . 1 17  *WK**2 . 

33650  BI-EXP (4 . *BX)-1 . 

33660  PVR-PVR+BI*(AC0S(Y)-Y*(l.-Y**2.)**0.5)*Y*0.0i 

33670  62  CONTINUE 

33680  PVR-PVR*16./PI 

33690  90  VR-VR*VX 

33700  VRI-1 .-EXP(-2 . *VR**0 . 5 ) 

33710  IF(RR.LT. 0.001)  GO  TO  91 

33720  WV-EXP ( 4 . *VR )- 1 . 


33730  GD-PVR/WV 

33740  VRI»VRI*GD**0.5 

33750  91  CONTINUE 

33760  RETURN 

33770  END 
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33780 

33790  C 

33800  C 

33810  C 

33820  C 

33830  C 

33840  C 

33850  C 

33860 

33870 

33880 

33890 

33900 

33910 

33920 

33930 

33940 

33950 

33960 

33970 

33980 

33990 


34000 

34010 

34020 

34030 

34040 

34050 

34060 

34070 

34080 

34090 

34100 

34110 

34120 

34130 

34140 

34150 

34160  C 

34170  C 

34180  C 

34190 

34200 

34210 


APPENDIX  C 


SUBROUTINE  SINTL(IV) 


THIS  SUBROUTINE  IS  TO  CALCULATE  THE  VARIANCE  OF  INTENSITY  OR 
POWER  DUE  TO  TURBULENCE  FOR  SPHERICAL  WAVE  SOURCE  OR  BEAM  WAVE 
WITH  PARTIALLY  COHERENT  SOURCE 

THE  CALCULATED  STANDARD  DEVIATION  CAN  BE  USED  TO  DEFINE 
HIGH  BOUND  AND  LOW  BOUND  OF  TRANSMITTANCE 


COMMON  / CARD 1 /  MODEL .IHAZE , ITYPE ,LEN , JP ,IM ,Ml ,M2 ,M3 ,ML .IEMISS ,K0 
1  , TBOUND , ISEASN , IVULCN , VIS 

COMMON  / CARD2/  HI ,H2 .ANGLE, RANGE, BETA, HMIN, RE, TT,RR,RLS 
COMMON  / CARD3/  VI , V2 ,DV ,AVM ,C0 ,CW,W( 15  )  ,E ( 15 ) ,CA,PI 
COMMON  / CNTRL/  LENST.KMAX ,M,IJ ,J1 ,J2 , JMIN,JEXTRA,IL,IKMAX,NLL,NP1 
1 , IFIND ,NL ,IKLO 

COMMON  /WANG/  K2 ,DSW(34 ) ,DZW(34 ) ,XWl (34 ) ,XW2 (34 ) 

COMMON  /VRAN/  VRI 
DIMENSION  WS(34) 

HW»H 1*1000.0 
WRANGE-RANGE*1000 .0 
WH2»H2* 1000.0 
GAA*0. 

TL0=*0. 


TLX=0. 

TLG-O. 

CN2*0.0 

WL-0.01/IV 

FR- (WL*WRANGE ) **0 . 5 

WK»IV*100.*2.*PI 

VK-WK**(7./6. ) 

DC-TT/RLS 

WK2«WK**2. 

WR2«WRANGE**2 . 

AM-WK/WRANGE 
SA-AM**2 . 

IF (RLS.LT. 0.0001)  RLS-0.0001 
RW«1./RLS**2. 

STT«TT**2 . 

IF (ITYPE .NE. 1 )  GO  TO  20 


VARIANCE  CALCULATION  FOR  HORIZONTAL  PATH 


CN2-4.2E-14*EXP(-HW/320.)/HW**0.667 
IF(Hl.GT.lOO.O)  CN2-0.0 
IF (HW.LE . 10.0 )  CN2-8.77E-15 


34220  SIGM-0.124*CN2*VK*WRANGE**1.833 

34230  RL0«(0.546*CN2*WK2*WRANGE)**1.2 

34240  RLX-0.425*CN2*WK**2.167*WRANGE**0.833 

34250  RLXS»0.114*RLX/0.425 

34260  GO  TO  80 

34270  20  IF (ANGLE. GT. 90.0)  GO  TO  37 

34280  C 

34290  C  VARIANCE  CALCULATIONFOR  UPWARD  PATH 

34300  C 

34310  DS-0.1 

34320  DO  34  I-J1.J2 

34330  WS(I)-DSW(I)/DZW(1) 

34340  IF (XWl(I).LE.25.0)  GO  TO  21 

34350  IF(XW1(I).LE.50.0)  GO  TO  22 

34360  DH-400.0 

34370  GO  TO  23 

34380  21  DH-20.0 

34390  GO  TO  23 

34400  22  DH-100.0 

34410  23  WIK«(XW2(I)-XW1(I))*1000.0/DH 

34420  IK-WIK 

34430  HW-XW1(I)*1000.0 

34440  DO  33  J-l.IK 

34450  HW-HW+DH 

34460  IF (HW.GE. 99600.0. AND. ITYPE.EQ. 3)  XW2 ( I )= 100000.0 

34470  IF(HW.GE. 99600.0)  GO  TO  34 

34480  CN2-4.2E-14*EXP(-HW/320. )/HW**0.667 

34490  IF(HW.LE.IO.O)  CN2-8.77E-15 

34500  DQ-WRANGE-DS 

34510  TLG*TLG+CN2  *DS  *  *0 . 8 3  3 *DH*WS ( I ) 

34520  TLO-TLOfCN2*DS**l .667*DH*WS (I ) 

34530  TLX-TLX+CN2*DQ**1.833*DH*WS(I)/DS**0.167 

34540  DS-DS+DH*WS ( I ) 

34550  33  CONTINUE 

34360  DQ-WRANGE-DS 

34570  IF(DQ.LE.O.)  DS-DS-DH*WS (I ) 

34580  IF(DQ.LE.O. )  DQ-DQ+DH*WS(I) 

34590  XY-(XW2(I)*1000.-HW)*WS(I) 

34600  IF(XY.LE.O.)  GO  TO  34 

34610  TLG-TLG+CN2*DS**0.833*XY 

34620  TLO-TLOK:N2*DS**1.667*XY 

34630  TLX-TLX+CN2  *DQ**1 . 833*XY/DS  **0 . 1 6  7 

34640  DS-DS+XY 

34650  34  CONTINUE 


K  ••:■ 

34660 

GO  TO  75 

34670 

37 

CONTINUE 

1 

34680  C 
34690  C 

VARIANCE  CALCULATION  FOR  DOWNWARD  PATH 

jnv 

34700  C 
34710 

DS-0. 1 

R£ 

34720 

DT-0 . 1 

34730 

Ll-Jl 

34740 

DO  60  L-l.NL 

34750 

WS (Li )*DSW(L1 )/DZW(Ll ) 

<W 

34760 

IF (XW1 (LI ) .LE .25 .0 )  GO  TO  38 

,  A  , 
•  ,  < 

34770 

34780 

IF (XW1 (LI ) -LE. 50.0 )  GO  TO  39 

DH-400.0 

>> 

34790 

GO  TO  40 

-\V 

34800 

38 

DH-20.00 

* 

34810 

GO  TO  40 

w 

34820 

39 

DH-100.00 

34830 

40 

WIK-  (XW1  (LI  )-XW2  (LI  )  )*1000 . 0 /DH 

34840 

IK-WIK 

34850 

HW-XW1(L1)*10Q0.0 

34860 

DO  57  J-l.IK 

m 

3487C 

HW-HW-DH 

i 

34880 

CN2-4 . 2E-14*EXP (-HW/320 . )/HW**0 . 667 

t® 

34890 

IF(HW.LE.IO.O)  CN2-8.77E-15 

f8S 

34900 

DQ-WRANGE-DS 

K' 

34910 

TLG-TLG+CN2*DS**0.833*DH*WS  (LI ) 

m 

34920 

TLO=TLO+CN2  *DS  ** 1 .667*DH*WS (LI ) 

34930 

TLX-TLX+CN2*DQ**1 .833*DH*WS  (LI  )/DS**0. 167 

34940 

34950 

DS-DS+DH*WS(L1) 

IF(K2.EQ.0)  GO  TO  57 

£>: 

34960 

IF (K2 .EQ . 1 .AND. WH2 .LE.HW )  GO  TO  57 

*  •  " 

34970 

TLG-TLG+CN2* (WRANGE-DT )**0 . 833*DH*WS (LI ) 

v:;:' 

34980 

TL0-TLO+CN2* (WRANGE-DT)**1 . 667*DH*WS (LI ) 

34990 

TLX-TLX+CN2*DT**1 ,833*DH*WS (LI )/ (WRANGE-DT)**0. 

35000 

DT-DT+DH*WS (LI ) 

K-'- 

35010 

57 

CONTINUE 

35020 

DQ-WRANGE-DS 

SB! 

35030 

IF (DQ .LT.O . )  DQ-DQ+DH*WS (LI ) 

ss 

35040 

IF(DQ.LT.O.)  DS-DS-DH*WS (LI ) 

SB 

35050 

XY-(HW-XW2 (LI ) *1000. )*WS (LI ) 

lx 

35060 

IF (XY . LE.O . )  GO  TO  61 

35070 

TLG-TLG+CN2  *DS  **0 . 8  3  3  *XY 

!«§ 

35080 

TLO-TLO«:N2*DS**1 .667*XY 

tgi 

35090 

TLX-TLX+CN2*DQ**1 ,833*XY/DS**0 .167 

v>'-> 


J-V*. 

fife 


A.-I'aV^a'a 


'-»•-*  ■ 


•.'  *-*~  •'■  •»'  .'-'.  AA -'■  -' ' -*■  -'■- u'uVj 


35100  DS-DS+XY 

35110  IF(K2.EQ.O)  GO  TO  58 

35120  IF(K2.EQ.1.AND.WH2.LE.HW)  GO  TO  58 

35130  TLG-TLG+CN2*(WRANG£-DT)**0.833*XY 

'35140  TLO-TLOfCN2*(WRANGE-DT)**1.667*XY 

35150  TLX-TLX+CN2*DT**1.833*XY/(WRANGE-DT)**0.167 

35160  DT-DT+XY 

35170  58  CONTINUE 

35180  Ll-Ll-1 

35190  IF (K2.EQ.0.AND.L1.LE.J2)  GO  TO  61 

35200  IF (LI .LE.JMIN.AND.K2 .EQ . 1 )  GO  TO  61 

35210  60  CONTINUE 

35220  61  CONTINUE 

35230  75  CONTINUE 

35240  SIGM«0.225*VK*TLG 

35250  RLO-1 .57rj*WK**2.4*TLO**l  .2/WRANGE**2 . 

35260  RLXS“0.235*WK*VK*TLX/WRANGE**1 .833 

35270  80  CONTINUE 

3528t»  F4-EXP(4.*SIGM) 

35290  C 

35300  C  CALCULATION  FOR  SPHERICAL  WAVES 

35310  C 


35320  IF(TT.GE. 0.001)  GO  TO  81 

35330  F5-1 .+RR**2 . *SIGM*RL0*8 . 

35340  SMP-F4 / F5  * (2 . +EXP ( -4 . *GAA ) -2 . *EXP ( -2 . *GAA )  )- 1 . 

35350  GO  TO  89 

35360  81  CONTINUE 

35370  C 

35380  C  CALCULATION  FOR  BEAM  WAVES 

35390  C 

35400  RDI-l./STT 

35410  F1«4.*RDI**2. 

35420  TM-2 . *RDI+2 . *RW 

35430  XM-0.5*RDI+0.25*RW 

35440  HX«RL0*SIGM*2 . 

35450  AML-AM-RLXS*4. 

35460  AML S* AML* *2. 

35470  BM-0.25*TMfSTT*SA/2.+4.*RL0 

35480  DM-0.5*RDI+HX*2. 

35490  EK-DM+0.5*RW 

35500  GM-EK+AMLS/(4.*DM) 

35510  EN-DM+AMLS/(4.*EK) 

35520  QM-DMf0.25*RW 

•35530  UM-SA/TMf0.25*TM+4.*RL0 
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35540  WM-GM-0.5*RW 

35550  YM-XM+SA/ (4 • *XM)+4 . *RL0 

35560  ZM-QM+AMLS/(4.*QM) 

35570  HA-SA/BM 

35580  HB-SA/UM 

35590  HC-SA/YM 

35600  F2-(TH/8.0+RLCH-STT*SA/4.  )**2. 

35610  FA-0 . 5  *STT / (BM*DM*GM ) 

35620  FB-0.5*STT/(BM*EK*£N) 

35630  FC-EXP(-4.*GAA)/ (TM*U M*DM*WM ) 

35640  FD-EXP(-2.*GAA)/ (2 .*XM*YM*QM*ZM) 

35650  F3-FA+FB+FC-FD 

35660  SMI-Fl*F2*F3*F4-l . 

35670  IF (SMI .LE .0 . )  SMI-0. 

35680  VRI-SMI**0.5 

35690  IF(RR.LT. 0.001)  GO  TO  90 

35700  SF1-HX*2 .-HX**2 . /DM+(AML*HX/DM4-AM+RLXS*2 . )**2 . / (GM*4 . ) 

35710  SF2-HX*2 . -HX**2 . /EK+(AML*HX/EK+AMfRLXS*2 .)**!./ (EN*4 . ) 

35720  SF3-HX*2 .-HX**2 . /DM+  <  AML  *HX / DM+AM+RLXS  *2 . )**2 . / (WM*4 . ) 

35730  SF4-HX*2 . -HX**2 . / QM+ (AML*HX/ QM+AMf  RLXS  *2 . )**2 . / (ZM*4 . ) 

35740  Q1-F1*F2*F4 

35750  Q2-STT*SA/ (1 .+DC**2 .+STT**2 .*SA+STT*4 . *RLO) 

35760  RR2-1./RR**2. 

35770  Q3-FA/ ( (SFl+RR2*0 . 5 )*(HA+RR2*2 . ) ) 

35780  Q4-FB/ ( (SF2+RR2*0 . 5 ) * (HA+RR2*2 . ) ) 

35790  Q5-FC/ ( (SF3+RR2*0 . 5 )*(HB+RR2*2 . ) ) 

35800  Q6-FD/ ( (SF4+RR2*0. 5 )*(HC+RR2*2 . ) ) 

35810  SMP«Ql*(Q2+RR2)**2.*(Q3+q4-K}5-Q6)-l. 

35820  89  IF (SMP.LE.O. )  SMP-O. 

35830  VRI-SMP**0.5 

35840  90  CONTINUE 

35850  RETURN 

35860  END 
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Plane  wave  sources;  point  receivers;  horizontal  paths. 


Fig.  2.  Plane  wave  sources;  point  receivers;  upward  paths. 
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Beam  wave  sources,  a  =  2cm;  point  receivers;  downward  paths 


»*<• 


zontal  pat 


Fig.  12.  Beam  wave  sources,  a  *  2cm;  receiver  radius  R  =  10cm;  downward  paths 
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Fig.  14.  Spherical  wave  sources;  point  receivers;  upward  paths 


Fig.  16.  Spherical  wave  sources;  receivers  radius  R  =  10cm;  horizontal  paths. 


Fig.  17.  Spherical  wave  sources;  receicers  radius  R  =  10cm; upward  paths. 
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MISSION 
of 

Rome  Air  Development  Center 

KADC  plan*  and  execute*  research,  development,  test  and 
4e£ec ted  acquisition  programs  in  support  of  Cotmand,  Control 
Communications  and  Intelligence  1C* I)  activities.  Technical 
and  engineering  support  within  areas  of  technical  competence 
as  provided  to  ESP  Program  0  a  ices  IPOs]  and  other  ESO 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur- 
veUCance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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